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Abstract For given nonnegative constants g, h, p, a with p 2 + a 2 = 1 and g + 
h > 0, we construct a diffusion process (Xi(-),X2(-)) with values in the plane and 
£3 infinitesimal generator 
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We compute the transition probabilities of this process, discuss its realization in terms 
of appropriate systems of stochastic differential equations, study its dynamics under a 
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driven by suitable continuous, orthogonal semimartingales M(-) and N(-) and with 
/(•) of bounded variation, which we study here in detail; and those of a one-di- 
mensional diffusion Y(-) with bang-bang drift dY(t) = —A sgn(Y(?))df +dW(t), 
7(0) = y driven by a standard Brownian motion W(-) . 

We also show that the planar diffusion (Xi (-),X2(-)) can be represented in terms 
of this process 7(-), its local time L Y (■) at the origin, and an independent standard 
Brownian motion Q(-), in a form which can be construed as a two-dimensional ana- 
logue of the stochastic equation satisfied by the so-called skew Brownian motion. 

Keywords Diffusion • local time • bang-bang drift • Levy characterization 
of Brownian motion • Tanaka formulae • weak and strong solutions • skew 
representation • skew Brownian motion • modified and perturbed Tanaka equations • 
time reversal 
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1 Introduction 

For given nonnegative constants g, h, p , a with p 2 + a 2 > and g + h > 0, and 
for a given vector (x\,X2) € K 2 , we shall consider the question of constructing a 
two-dimensional diffusion process (Xi(-),X2(-)) with dynamics 

dXi (f) = (gl{xi(0<x 2 (0} ~ hl {xdt)>x 2 (t)}) dt + (P 1 {x l (t)>x 2 (t)} + ^1{X!(0<x 2 (i)}) dfii (f ) , 

(1.1) 

dX 2 (0 = (gl{xi(t)>x 2 (t)} - hl {Xi(t)<x 2 (i)}) dt + (P 1 {x l (t)<x 2 (t)} + ^{Xi(t)>x z (t)}) dB 2 {t) , 

(1.2) 

initial condition (X\ (0),X2(0)) = (x\ ,x 2 ), and B\ (•), B2G) two independent, standard 
Brownian motions. For simplicity, we shall use throughout the normalization 

p 2 + a 2 = l (1.3) 

and refer to the case pa = as "degenerate". 

Speaking informally and a bit imprecisely for the moment about the system of 
(1.1)-(1.2), imagine you run two Brownian-like particles on the real line. At any given 
time, you assign positive drift g and diffusion a to the laggard; and you assign 
negative drift —h and diffusion p to the leader. What is the probabilistic structure 
of the resulting two-dimensional diffusion process? Can it be realized as the solution 
of a system of stochastic differential equations other than (1.1), (1.2)? What are its 
transition probabilities? How does it look like, when time is reversed? 

It has been known for some time now, at least for the non-degenerate case (cf. 
Stroock & Varadhan (1979), pages 193-194; Bass & Pardoux (1987); or 
Krylov (2004), page 45) that a unique probability measure fl can be constructed 
on the canonical filtered measurable space (2XJ,5f), G = {^(f)}o<r<~ of continuous 
functions tn : [0,°°) — ?► M 2 endowed with the topology of uniform convergence on 
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compact intervals, such that the process /(ro(f)) — / jSf/(tn(i))ds, < t < °° is a 
local martingale under /i for every function / : R 2 — >• K of class 'iS 1 (here J?f is the 
second-order partial differential operator in (0.1), the infinitesimal generator of the 
resulting diffusion). Our goal in this paper is to describe this probability measure ju 
as explicitly as possible; to study its behavior under time-reversal; and to understand 
the solvability of systems of stochastic differential equations, such as (1.1), (1.2) 
above, which correspond to this martingale problem and help realize its solution. 

We shall show in sections 2, 3 and 4 that the system of stochastic differential 
equations (1.1), (1.2) has a solution which is unique in the sense of the probability 
distribution (thus the above martingale problem is indeed well-posed). This solution 
is shown to be strong in section 5; it is characterized in terms of a one-dimensional 
diffusion process Y(-), which has "bang-bang" drift with intensity A = g + h > 
and is driven by yet another standard Brownian motion process W(-) , namely 



Here and in what follows we shall use the convention for the signum function 



The one-dimensional diffusion of (1.4) was studied in some detail by KARATZAS 
& Shreve (1984), who found its transition probabilities as in equations (6.3), (6.4) 
below. We shall use this analysis to compute, in section 6, the transition probabilities 
of the two-dimensional process (X\(-),X2(-)) . As in that earlier paper, a crucial role 
will be played here again by the local time 



accumulated at the origin during the interval [0,f] by the diffusion Y(-) of (1.4). In 
terms of this diffusion, its local time (1.6), and an independent standard Brownian 
motion Q(-), the unique-in-distribution weak solution of the system in (1.1), (1.2) 
will be shown to admit the skew representation of (2.22), (2.23) below. 

The diffusion process of (1.4) is also instrumental in a time-reversal analysis we 
carry out in section 7, where the dynamics of time-reversed versions of the pro- 
cesses Xi(-), Xz(-) are derived in the spirit of Haussmann & Pardoux (1986) 
and MEYER (1994). We were quite surprised, at first, that these reverse-time dynam- 
ics should feature terms involving singularly continuous components such as local 
times, as indeed they do. To the best of our knowledge, this is the first instance where 
such a structure is observed in the context of a "purely forward" system of stochastic 
differential equations such as (1.1)-(1.2) that does not involve reflection; see Remark 
7.3 in this regard. We also study the forward and backward dynamics of the ranks 
R\(-) = msK(Xi(-),X 2 (-)), Z?2(-) = nun(Xi('),Z 2 (-)) in this two-dimensional diffu- 
sion, in sections 4 and 7, respectively. 

The planar diffusion process with infinitesimal generator (0.1) has local covari- 
ance matrix 



dY(t) = -Asgn(F(f))df+dW(f). 



(1.4) 



sgn(y) :=l( ,=o)(y)-l(-=o,o]0), yeK. 



(1.5) 





(1.7) 



4 



E. Robert Fernholz et al. 



There is a continuum of real square roots for this matrix, of the form E(xi,X2) — 
ifa >.v 2 } + l{*i <* 2 } wim 

£ / pcos<p — psin<p\ / (TcosiJ — csin# 

+ ' yeasinp ecrcosp y 5 ~ yp5siniJ p5cos# 

parametrized by e = ±l, 5 = ±1, < <p,# <1%. All such configurations lead 
to systems of stochastic differential equations that admit a (unique in distribution) 
weak solution. We show in subsection 5.1 that those solutions that correspond to 
configurations with 

(ff 2 £-p 2 5)sin(#-<p)+pc7(l+£(5)cos(#-<p) = -1 (1.9) 

are not strong; see the system (2.14), (2.15) for an example. Whereas all other con- 
figurations lead to strongly solvable systems; one such system appears in (1.1), (1.2), 
and another one in (2.12), (2.13). 

Crucial in this analysis of strength and weakness is the following recent result by 
Prokaj (201 1) on the pathwise uniqueness of the "perturbed Tanaka equation" 

r(f) =y + J sgn(Y(s))6M(s)+N(t), 0<f<°o. (1.10) 

Theorem 1.1 (PROKAJ (201 1)). Suppose that M(-) , N(-) are continuous local mar- 
tingales with M(0) = A^(0) = and quadratic and cross-variations that satisfy the 
conditions of orthogonality and domination 

(M,N)(t) = 0, (M)(t) = [ q(s)d(N)(s); 0<f<°o, (1.11) 

Jo 

respectively, for some progressively measurable process q(-) with values in a com- 
pact interval [0, c] . Under these assumptions, pathwise uniqueness holds for the per- 
turbed Tanaka equation (1.10). 

In section 8 we shall use the local time techniques introduced by PERKINS (1982) 
and further developed by Le Gall (1983), to provide a simple proof of a consider- 
ably more general result of this type, Theorem 8.1 (see also Proposition 8. 1), in which 
the signum function is replaced in (1.10) by an arbitrary function of finite variation, 
and M(-) , N(-) by continuous semimartingales such that (1.1 1) is satisfied. 

Multidimensional processes of the type (1.1), (1.2) were introduced by FERN- 
HOLZ (2002), and their ergodic behavior was studied by BANNER et al (2005), PAL 
& Pitman (2008), Ichiba et al (2011), among others. Here we focus on the two- 
dimensional case, and concentrate on the precise probabilistic structure of the result- 
ing diffusions governed by stochastic differential equations such as (1.1), (1.2). 
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2 Analysis 

Let us assume that the system of stochastic differential equations (1.1), (1.2) has a 
weak solution: To wit, that there exists a filtered probability space (£2,$,W), F = 
{5 r (0}o<r<~ and on it two pairs (Bi(-),B 2 (-)) an d (-^i(')>^2(*)) of continuous, 
F— adapted processes, such that Bi(-) and B 2 (-) are independent standard Brow- 
nian motions and (1.1), (1.2), ^i(O) = X\ and ^(0) = X2 hold. We shall fix the 
nonnegative constants g, h, p, a with g + h > 0, and impose the normalization 
(1-3). 

We shall assume throughout the paper, and without further mention, that the fil- 
trations we are dealing with are in their right-continuous versions and have been 
augmented by sets of P— measure zero. We shall also use the convention 

F S = {3 S (0}<xr<~, S s (r) := a(S(s),0<s<t) 

for the P— augmentation of the filtration generated by a given process E : [Q,°°) x 
£2 — > W 1 with values in some Euclidean space and RCLL paths. 
With such a setup, and with the notation 

X = g + h, V = g — h, y = x\—X2, z = xi+x2, (2.1) 

we note that the difference 

Y(t):=X 1 (t)-X 2 {t), 0<f<°o (2.2) 

satisfies the integral version 

Y(t) = y-X f sgn(Y(s))ds + W(t), < t < °o (2.3) 

of the equation (1.4). The equation (2.3) is driven by the process 

W(t) := pWi(t) + aW 2 {t), 0<t<oo. (2.4) 

where we have set 

Wi(t) := f l {yW>0} dSi(*)-^ l {Y{s)m dB 2 (s), (2.5) 

W 2 (t) := f l^^ojdfliW- jT l {Y{s)>0} dB 2 (s). (2.6) 

It is also seen from (1.1) and (1.2) that the sum of the two processes Xi(-), X 2 (-) is 
of the form 

Xi(t)+X 2 (t) =z + vt + V(t), V{t) := pV x {t) + aV 2 {t), (2.7) 
where, by analogy with (2.5), (2.6) we have set 

V X (f ) := jT 1 { k(. s)> o} dB, (*) + £ 1{y(s)<0} dB 2 (s) , (2.8) 

V 2 (t) := f 1 { Y( S) < 0} dfli (*) + f 1 { y( s)>0} dB 2 (s) . (2.9) 
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• The processes Wi(-), W 2 (-), Vi(-), V 2 (-) are continuous F— martingales with 
quadratic variations (Wi)(t) = (W 2 )(t) = (Vi)(f) = (V 2 )(t) = t, i.e., Brownian mo- 
tions by the P. Levy theorem (e.g., Karatzas & Shreve (1991), p. 157); also 
(Vi,V 2 )(-) = (W h W 2 )(-) = (Vt,W 2 )(-) = (Wi,V 2 )(-) =0. 

The pairs (Wi(-),W 2 {-)) and (V^-), V 2 (-)), as well as (Vi(-),W 2 {-)) and (Wi (•),%(•)), 
are thus two-dimensional Brownian motions, by the P. Levy theorem. In conjunction 
with (1.3) this implies, in particular, that W(-) and V(-) in (2.4), (2.7) are standard, 
one-dimensional Brownian motions. 

• On the other hand, we see from (2.5), (2.8) and with the notation of (1.5) the 
intertwinements 

Vi(t) = J sgn(Y(s))dWi(s), Wi(t) = J sgn(r(*))dVi(*); (2.10) 
and from (2.6), (2.9) we obtain the intertwinements 

V 2 (t) = -J\gn(Y(s))dW 2 {s), W 2 (t) = -J\gn(Y(s))dV 2 (s). (2.11) 

• The equation (2.3) admits a pathwise unique, strong solution; in particular, unique- 
ness in the sense of the probability distribution holds as well (cf. Proposition 5.5.17, 
page 341 and the YAMADA-WATANABE Proposition 5.3.20, page 309 in KARATZAS 
& Shreve (1991)), and we have the identity F F = F w . 



2.1 Two Auxiliary Systems 

The relations X x (t) +X 2 (t) = x\ +x 2 + (g-h)t + pVi (t) +oV 2 {t) from (2.7), and 
Xx(t)-X 2 (t) = 

x\ -x 2 + (g+h) J (l{ Xl (s)<x 2 {s)} ~ 1 {x I (,)>x 2 ( s )}) ds + pWi(t) + aW 2 {t) 

from (1.4), (2.2), (2.4) lead, in conjunction with the intertwinements of (2.10) and 
(2.11), to the system of stochastic integral equations 



M*)=xi + f o «<jr 2 (*)} -fH{ Xl (s)>x 2 { t )}) 

+ pJ Q l{Ms)>Ms)} dw i( s ) + a J l ms)<x 2 (s)}dW 2 (s) 

-pf Q hM s) <x 2{ s) } m (s) cj jf' l {m>m} dW 2 (s) (2.13) 
driven by the planar Brownian motion (Wi (•))^2(')j °f (2-5), (2.6). 



ds 

(2.12) 
ds 
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Furthermore, we observe that the intertwinements of (2.10), (2.11) allow us to 
recast this system as driven by the planar Brownian motion (Vi(-),\^(«)) of (2.8), 
(2.9), namely 

*i(0 =Xl+ j Q (^ 1 WM<^W}- /z %iW>x 2 W})^ 

+ pJ Q hx^x^ydV^s) + oJ Q l {Xl (s)<x 2 (s)}dV2(s) , (2.14) 
X 2 {t) = X 2 + J Q (8i{x^)>x 2 (s)}- h^xas)^)}) ds 

+ P J^{xds)<x 2 (s)}dVi(s) + a j\ { x l{s> x 2{s ) } dV 2 (s) . (2.15) 

It is quite clear, though perhaps worth noting, that the systems of stochastic equations 

(1.1)-(1.2), as well as (2.12)-(2.13) and (2.14)-(2.15), give rise and correspond to the 
same martingale problem - namely, the one with infinitesimal generator (0.1). 

2.2 Skew and Integral Representations 

By analogy with (2.4), (2.7) let us introduce the standard Brownian motions 

w\-):=pWi(-)-aW 2 (-), V b (-):=PVi(-)-ffV 2 (-) (2.16) 
and note the new intertwinements 

V(t) = J sgn (Y(s)) dW''(s) , V\t) = sgn (Y(s)) dW(s) . (2.17) 

It will be convenient to cast the Brownian motion W (•) of (2.16) in the decomposi- 
tion 

w\-) = yW(-) + 8u\-), where y:=p 2 -a 2 , 8:=^/l-y 2 = 2pa , (2.18) 

in terms of the independent Brownian motions 

U b {-):=aWi{-)~pW 2 {-) and W(-) = pWi{-) + aW 2 {-) (2.19) 

as in (2.4). With this setup, the Brownian motion V (•) defined in (2.7) takes the form 
V(t) = Jo' sgn (y(s)) [ydW(s) + 8dU b {s)} , equivalently 



V(t) = yjf sgn (70)) [dF(i)+Asgn(r(»)ds] +5 Q(t) 

= y(\Y(t)\-\y\+kt -2L Y (t))+8Q(t), 0<f<oo (2.20) 



thanks to (2.17), (2.18), (1.4) and the Tanaka formulas (e.g., Karatzas & Shreve 
(1991), page 220). Here 

Q{t):= J sgn (Y(s))du\s) = oVi{t)+pV 2 {t), 0<r<- (2.21) 
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is yet another Brownian motion with (Q,W) = (U ,W) = 0. Thus, Q(-) is inde- 
pendent of the Brownian motion W(-) and of the diffusion process Y(-) that W(-) 
engenders via (2.3). 

• We can combine now the last expression in (2.20) with X\(t) — X 2 (t) = Y(t) of 
(2.2) and with X\{t) +X 2 {t) = x\ + x 2 + Vt + V(t) of (2.7), and arrive at the skew 
representations 

X l {t)^x l +iit+p 2 {Y + {t)-y+)~a 2 {Y-{t)~y-)~yL Y {t)+paQ{t), (2.22) 

X 2 (t) =x 2 +llt- a 2 (Y + {t)-y+) +p 2 ( Y -(t)-y-) - yL Y ' (t) + poQ(t) (2.23) 

for the components of the two-dimensional diffusion (Xi(-),X 2 (-)) of the system 
(1.1), (1.2); we have set 

li:= l -{v + Xy) = gp 2 -ha 2 . (2.24) 

These formulas involve the positive and negative parts of the current value of the 
one-dimensional diffusion process Y(-) in (1.4), the current value of its local time 
L Y (•) at the origin, and the current value of the independent Brownian motion Q(-) 
of (2.21). Two cases stand out. 

(A) In the Equal Variance (Isotropic) case p = a = \ /v2 , the local times disappear 
from the expressions (2.22), (2.23), which then take the very simple form 

Xi (t) = xi + i (vt + Y(t) -y + Q(t)) , X 2 (t) =x 2 + \ (Vt - Y(t) +y + Q(t)) . 

(2.25) 

(B) In the Degenerate case pa = the independent Brownian motion Q(-) disap- 
pears from these expressions; for instance, with (7 = and p = 1 , they become 

Xj.it) =x l -y + +gt+Y + (t)-L Y (t), X 2 (t) =x 2 -y-+gt+Y-(t)~L Y (t). 

(2.26) 

2.3 Uniqueness in Distribution 

The analysis of this section shows that, given any weak solution of the system of 
stochastic equations (1.1), (1.2), its vector process (Xi(-),X 2 (-)) can be cast in the 
form (2.22)-(2.23). Here the diffusion process Y(-) is the pathwise unique, strong 
solution of the stochastic integral equation 

Y{t)=y-X j sgn(Y(s))ds+W(t), < t < °o 

as in (2.3) and with the notation of (1.5), driven by the Brownian motion W(-) of 
(2.4); whereas the Brownian motion Q(-) is independent of W(-), thus also of Y(-). 
In other words, the joint distribution of the pair {Y (■),()(■)) is determined uniquely 
- and thus, from (2.22)-(2.23), so is the joint distribution of the vector process 

(*l(0>*2(-)). 

To put it a bit more succinctly: uniqueness in distribution holds for the system 
of equations (1.1), (1.2), as well as for the systems of equations (2.12), (2.13) and 
(2.14), (2.15). 
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3 Synthesis 

Let us begin now to reverse the steps of the preceding analysis. We start with a filtered 
probability space (Q,$,F), F = {5tf)}o<r<°° rich enough to support two indepen- 
dent, standard Brownian motion Wi(-) and W 2 (-) and without sacrificing generality, 
we shall assume F = F( w i' W 2\ i.e., that the filtration is generated by this planar Brow- 
nian motion. 

With given nonnegative constants g, h, p, O that satisfy (1.3) and g + h > 0, 
with a given vector (xi,x 2 ) £ K 2 , and with the notation of (2.1), we construct the 
pairs of independent Brownian motions 

W(-) : = pW 1 (-) + aW 2 (-), U b {-) : = aWi(-) - pW 2 {-) (3.1) 

and 

[/(•) := dWi(-) +pW 2 {-), W b {-) : = pWi(-) - aW 2 {-) (3.2) 

in accordance with (2.19), (2.16). Clearly, F^^ eF^) = F^* . 

We construct also the pathwise unique, strong solution Y(-) of the stochastic 
equation (2.3) driven by the Brownian motion W(-) in (3. 1). This is a strong MARKOV 
and FELLER process, whose transition probabilities can be computed explicitly; see 
(6.3)-(6.5) below. 

With the process Y(-) thus in place, we introduce the continuous, F— adapted 
processes 

Vi(r) - J sga(Y(s))dWi{s), V 2 {t) - - J sgn(Y(s))dW 2 (s) (3.3) 

in accordance with (2. 10) and (2.1 1). These are martingales with (Vi)(f ) = (V2) (0 = 
t and (Vi,V2)(f) =0 for all < t < °o, thus independent Brownian motions in their 
own right. We construct from them, and by analogy with (3.1) and (3.2), two addi- 
tional pairs of independent, standard Brownian motions, namely 

V(-): = pVi(-) + aV 2 {-), Q"{-) ■ = aVi(-) - pV 2 (-) (3.4) 

and 

g(-) := aVi(-)+pV 2 (-), v\-): = pV x (-)-aV 2 {-). (3.5) 

We note the intertwinements (2.17), (2.21) and Q\-) = J Q sgn (Y(t )) dU(t ) , as well 
as the filtration identities F^'^ =F( v -2 b ) = F^ b ) . 

Finally, we introduce the continuous, F— adapted processes 

Zi(t) : = xi + J (gl {Y{s)m - hl {Yis)>0} ) ds + Mi(t) (3.6) 
X 2 (t): = x 2 +f (gl {Y{s)>0} -hl [Y{s)m )ds + M 2 (t) (3.7) 
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for < t < °o , where we set 

Mi (t):= jT(p 1{f(.v)>o } dWi (*) + a 1 {F( ,)< 0} m(sj) , (3-8) 
M 2 (f) : = £ (-pl Ws) < 0} dW 1 ( i )-(Jl {F(i . )> o } dW 2 (5)) , (3.9) 

by analogy with the equations of (2.12), (2.13). We have for these processes 

X 1 (t)-X 2 (t) = Y(t), X 1 (t)+X 2 (t) = Xl +x 2 + Vt + V(t) (3.10) 

in accordance with (2.2), (2.7), and note that M\(-), M 2 {-) are continuous F— mar- 
tingales, with (Mi,M 2 ) (•) = and quadratic variations 

(Mi){t) = J (p 2 l{y ( ,)>0} +CT 2 l Wj) < })ds , 
(M 2 )(t) = J (p 2 l {Y[s)m + a 2 l {Y{s)>0} )ds. 

There exist then independent Brownian motions Bi(-), B 2 (-) on our filtered proba- 
bility space (Q,$,F), F = {${t)}o< t< oo, so the continuous martingales of (3.8), (3.9) 
are cast in their Doob representations 

M\{t) = f (pl {F(s)>0} + (Tl Wj) <o})dBi(i), 

Jo t (3.11) 

M 2 {t) = (pl{F( s )<0} + <Tl{y( 4 .)>0})d£ 2 (i) 

for < f < °° ; for instance, we can take the Brownian motions 

2?i(f)= ^(l { y W> o } dWi(*) + l {yw < 0} dW 2 (5)), (3.12) 
ft(0 - - jT (l { y W < }dWi(5)+l { K W> o } dW 2 W) (3.13) 

that one gets by disentangling (Bi(-),B 2 (-)) from (Wi (•), W 2 (-)) in (2.5), (2.6). 

3.1 Taking Stock 

To recapitulate: we have constructed a weak solution for the system of stochastic 
differential equations (1.1), (1.2), as is seen clearly from (3.6), (3.7) and (3.11); and 
as we stressed in subsection 2.3, this solution is unique in distribution. We remarked 
in the Introduction that this is in accordance with general results of STROOCK & 
Varadhan (1979), pages 193-194; BASS & Pardoux (1987); or Krylov (2004), 
page 45. 
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It develops from (3.6)-(3.9) that we have also constructed a weak solution to the 
system (2.12), (2.13), once again unique in distribution. On the other hand, we can 
express the martingales of (3.8), (3.9) as 



in terms of the independent, standard Brownian motions of (3.3). Back into (3.6) and 
(3.7), these expressions show that we have also constructed a weak solution for the 
system of stochastic equations (2.14), (2.15), once again unique in the sense of the 
probability distribution. 

In the next two sections we shall discuss in detail properties of strength/weakness 
for the solutions to these systems of equations, namely, (2.14)-(2.15), (2.12)-(2.13) 
and (1.1)-(1.2). For the moment, let us remark that the Brownian motion V(-) = 
pVi(-) + aV 2 {-) determines the sum of Xi(-),X 2 {-) as in (2.7) or (3.10); and that 
the Brownian motion W(-) = pWi(-) + aW 2 (-) determines the difference Y(-) — 
Xi (•) — X 2 (-) via the solution of the stochastic equation (2.3). 

3.2 A Two-Dimensional Analogue of the Skew Brownian Motion 
The equations of (2.22)-(2.23) can be cast in the form 



namely, as a system of stochastic differential equations involving the local time L Xl ~ Xl (-) 
at the origin of the difference X\ (•) — X 2 (-) , We have constructed a solution to this 
system, subject to the requirement that the difference X\ (■) — X 2 (-) be independent of 
the driving Brownian motion Q(-). 

The system (3.14), (3.15) can be construed as a two-dimensional analogue of the 
stochastic equation derived by HARRISON & Shepp (1981) for the WALSH (1978) 
skew Brownian motion; see the recent survey by Lejay (2006), and the references 
cited there, for the various constructions and properties of this process. 

4 Ranks and Filtrations 

Let us look now at the ranked versions 





+ p 2 ((Xi(t)~X 2 (t))--(. 




a 2 )L x ^ x Ht)+paQ{t), 



max(Xi(-),X 2 (-)), 



min(Xi(-),X 2 (-)) 



(4.1) 
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of the components of the diffusion process constructed in the previous section. We 
have R l (t)+R 2 {t) =X l (t)+X 2 (t) =xi +x 2 + Vt + V(t) from (3.10), as well as 

R!{t)-R 2 (t) = |Xi(r)-X 2 (0| = \Y(t)\ = \y\+ J\gn(Y{s))dY{s) + 2L Y (t) 
= \xi-x 2 \ -It + V b {t) + 2L Y (t), 0<f<°o 

(4.2) 

from the Tanaka formulas and (2.3), (2.17); and from the theory of the SKO- 
rokhod reflection problem (e.g., Karatzas & Shreve (1991), page 210), we 
note 

2L Y (t) = max (- (\y\+v\s) - ls)^j , 0<f<°°. (4.3) 
We also note from (4.2), (4.3) and (1.6) the filtration relations 

d V \t)=$W(t)C.$ Y (t), <*<<*>; (4.4) 

the strict inclusion comes from the fact that Y(-) has a zero on [0,f] with posi- 
tive probability, so that the random variable sgn(F(f)) is not ^l F l (f)— measurable. 
Alternatively, one can note from (6.2) below that both P(F(f) > 0|S r l r l(f)) and 
P (Y (f ) <0\d m (f)) are non-trivial. 

With r\ = max(xi,X2), r 2 = min(xi ,x 2 ) we deduce from these equations and 
(3.4), (3.5) the dynamics for the ranks 

Ri(t) = n -ht + pVi{t)+L Y (t), 0<f<°o (4.5) 
R 2 (t) = r 2 +gt + oV 2 (t)-L Y (t), 0<r<oc. (4.6) 

Equations (4.5), (4.6) identify the processes Vi(-) and V 2 ( ) of (3.3) as the inde- 
pendent Brownian motions associated with individual ranks; such an interpretation 
is also possible using the equations of (2.14), (2.15). On the other hand, the inde- 
pendent, standard Brownian motions B\(-), B 2 (-) of (3.12) and (3.13), are those 
associated with the "names" (indices, identities) of the individual particles. 



4.1 The Degenerate Case 

We embark now on a detailed study of the solutions to the systems of equations 
(1.1), (1.2), as well as of those of (2.12), (2.13) and (2.14), (2.15). In this effort it is 
instructive to focus first on the degenerate case p a = , which has some interesting 
features of its own. 

Proposition 4.1. In the degenerate case with (7 = 0, thus p = 1 in light of (1.3), we 
have the relations 

for every < t < °°, where the inclusion is strict. 
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Proof. With p = 1 , a = we have Vi(-) = V(-) = V\-) in (3.4), (3.5); the first of 
the claims in (4.7) follows from (4.3) and (4.5)-(4.6), which read now 

R l {t) = r l -ht + V{t)+L Y (t), R 2 (t)=r 2 +gt~L Y (t), 0<*<°°. (4.8) 

The second claim is immediate from (3.10). 

The third (strict inclusion) is fairly obvious from (2.26). Indeed, the equations of 
(2.26) show 

X l {t)+X 2 {t) =z-\y\ +2gt + \Y(t)\ -2L Y {t) , 

thus $ Xl+X2 (t) C £l y l(?) C $ Y (t) for < t < oo; and, as we remarked in (4.4), the 
second inclusion is strict for < t < °° . 

The fourth and fifth relations in (4.7) follow from the equations in (2.26) and 
X\(-) —X 2 (-) = Y(-) ; whereas the sixth relation is a consequence of the strong solv- 
ability of the stochastic equation (2.3). □ 

The strictness of the inclusion in ^( fi '*'(f) = $ v (t) C $ w (t) = £( x i*)( f ) f or 
< f < oo reflects the fact that some information is inevitably lost when one passes 
from the "names" (X 1 (-),X 2 (-)) to the "ranks" (R\ (-),R 2 (-)) . In the case p = 1 , 
a = that we are studying here, the equations (2.12), (2.13) read 

hi-{ Xl ( s )>x 2 {s)}) d* + jf l {Xlis)>X2is)} dW(s) , 

(4.9) 

hl {x t (s)<x 2 U)}) ds - I 1 {X l (s)<X 2 (s)}dW(s) 

(4.10) 

with W(-) = Wi(-) in accordance with (2.4). The construction (synthesis) of section 
3 shows that the solution constructed there is strong, that is, adapted to the filtra- 
tion F w generated by the Brownian motion W(-) that drives this system. Repeating 
the analysis of section 2, we see that any solution (X\(-),X 2 (-)) of the system (4.9), 
(4.10) can be cast in the form (2.26), where Y(-) is the pathwise unique, strong so- 
lution of the diffusion equation (2.3) with F F = F w and L Y (■) is the local time in 
(1.6). In particular, any such solution of the system (4.9), (4.10) is a non-anticipative 
functional (given as in (2.26)) of its driving Brownian motion W(-) . Thus, pathwise 
uniqueness holds for the system of equations (4.9), (4.10). 

On the other hand, and always in the case p = 1 , a = , we can cast the equations 
(2.14), (2.15) as 

Xi(t) =*i+J o Ul{x,M<x 2 (.v)} -hl {Xl (s)>x 2 (s)})ds + j o l {Xl{s)>X2m dV(s), 

(4.11) 




(4.12) 

now driven by the Brownian motion V(-) = V\ (•) in accordance with (3.4). We know 
from Proposition 4. 1 that this system does not admit a strong solution. 



Xl(t) = X l+j Q Ul{X l( .v)<X 2 (.v)} 

X 2 (t) =*2 + J Q (sl{x,(.v)>x 2 ( s )} 



14 



E. Robert Fernholz et al. 



Collecting the results of this section and of the preceding one, we see that we have 
established the following theorem; this can be seen as a two-dimensional analogue of 
Barlow (1988). 

Theorem 4.1. The system of stochastic differential equations (1.1), ( 1.2) with p a = 
has a weak solution which is unique in the sense of the probability distribution. The 
same is true for each of the systems of equations (4.9), (4.10) and (4.11), (4.12). 

On the other hand, the system of stochastic differential equations (4.9), (4.10) 
admits a strong solution, which is therefore pathwise unique; whereas the system 
(4.11), (4.12) does not admit a strong solution. 

We have deduced here pathwise uniqueness from the "obverse Yamada & Watan- 
abe" results of Engelbert (1991) and Chernyi (2001). 



The following Figure is a Simulation of the Diffusion with p = 1 , a = and g = h = 1 . 




i 1 1 1 1 1 1 1 1 i r 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Fig. 1 Simulated processes; Black = X[ (■) , Red = Xj(-) . 



4.2 The Non-Degenerate Case 

Let us move now on to the non-degenerate case p a > . In contrast to the situa- 
tion we encountered in Proposition 4.1, now both pairs (Xi(-),X2(-)) (of positions 
by "name") and (/?i(-),i?2(0) (°f positions by "rank") generate two-dimensional 
Brownian filtrations, as our next result shows. 

Proposition 4.2. In the non-degenerate case p O > 0, we have for every < t < °° 
the filtration relations 

3(vi.v«(,) = $^> R i)(t) = d il¥lv) (t) - Z^\t) c 5 »*)( f ) = 

= 3 {Y ' v) (t) = d {w ^Ht) = d {¥ ' ub) (t) = s (w - e) (0 = z im (t), (4.i3) 
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where the inclusion is strict. 

Proof. We have clearly ^ (t) C £ yb (f ) C #W ^ (t) by virtue of (4.3), (3.4); back 
into the equations of (4.5), (4.6), this implies ^'^(t) C 3 iv ^ v i\t). For the re- 
verse inclusion, we note that S" L> '(f) C gl F l(f) C ^ (Rl ^(t) holds, thanks to (1.6) 
and \Y(-) \ ^Ri(-) -R 2 {-). Back into (4.5), (4.6) this gives & v ^{t) C S^'^O), 
and the first equality in (4. 13) is proved. 

In conjunction with (2.7) and (2.20), these considerations also give $( Rl,R2 \t) = 
g(\Y\,x l+ x 2 )^ = $(\Y\,V)(t) = $(\ Y \,Q)( t ) f or every < f < ~> , justifying the second 
and third equalities in (4.13). 

For the fourth equality in (4. 1 3), we have noted already the inclusion $( Xl ' X2 \t) C 
§( Wl ' Wl \t) from the construction of (3.6)-(3.9). In order to argue the reverse inclu- 
sion, let us note that (3.8), (3.9) imply 



On the strength of (3.6) and (3.7), the F^ 1 ^ -martingales Mj(-) and M 2 (-) are also 
martingales of the filtration F' Zl ' X2 ' generated by the state process (X\ (■), X2(-)) , so 
we deduce fW- 1 ^ C F**' and the fourth equality in (4.13) follows. 

For the fifth equality we note that a {X\(t) ,X 2 {t)) = o(Y(t ),V (f )) holds for ev- 
ery f > 0, courtesy of (3.10), and gives ^^(f) = 3 {Yy \t). The sixth equal- 
ity in (4.13) is a consequence of (2.19), which gives actually the stronger statement 
a(Wi (t),W 2 (t)) = <7(W(t),U b {t)) for every f >0. 

The seventh and eighth equalities in (4.13) are straightforward consequences of 
F w = F F and of the equation (2.21), along with its "twin" U b (•) = / sgn(F(f)) dg(f ) . 
The ninth equality is now obvious. 

Finally, the strictness of the inclusion in (4. 13) follows now from the strictness of 
the inclusion in 



itself a consequence of the strictness of the inclusion in (4.4) and of the independence 



Arguing as in the previous subsection and using Proposition 4.2, we obtain the 
following analogue of Theorem 4. 1 for the non-degenerate case. 

Theorem 4.2. The system of stochastic differential equations (1.1), ( 1.2) with p o > 
has a weak solution which is unique in the sense of the probability distribution. The 
same is true for each of the systems of equations (2.12), (2.13) and (2.14), (2.15). 

On the other hand, the system of stochastic differential equations (2.12), (2.13) 
admits a strong solution, which is therefore pathwise unique; whereas the system 
(2.14), (2.15) admits no strong solution. 




^*2) {t) = %{\n,Q)(t) c d (r,Q)(t) =$(*i.*2)(,) j <f 



of the processes Y(-), Q( ) . 



□ 
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The existence and uniqueness-in-distribution of a weak solution to (1.1), (1.2) 
in the non-degenerate case p a > follow also from the general, multidimensional 
results of BASS & Pardoux (1987). The approach here is more direct and con- 
crete, capitalizing on the special two-dimensional nature of our setting; the Bass- 
Pardoux results cannot, however, be applied to the degenerate case p a = 0, so it 
does not seem possible to obtain even the first part of Theorem 4. 1 from them. 

We shall see in Theorem 5.1 of the next section that the system of equations (1.1), 
(1.2) admits a strong solution; this too goes beyond the Bass & Pardoux (1987) 
results. As far as we know it is an open issue, whether such strength might obtain in 
their setting as well (to wit, where Mr , d > 2 is partitioned into disjoint polyhedral 
chambers, diffusion characteristics are constant in each chamber, and strong non- 
degeneracy prevails). 

4.3 Discussion 

Theorem 4.2 shows that the system of (2.14), (2.15) with p a > can be thought of 
as a genuinely two-dimensional TANAKA example: a system of stochastic differential 
equations that admits a weak solution which is unique in the sense of the probability 
distribution, but no strong solution. 

Theorem 4.2 might usefully be compared also with the strictly one-dimensional 
results by BARLOW (1988) and NAKAO (1972), as well as with Theorem 2 in Vere- 
tennikov (1982); see also Barlow (1982). Barlow (1988) shows, in particular, 
that with any given constants x E M , p > and a > , the one-dimensional stochas- 
tic differential equation 

X(t) =x + J z(X(s))6B(s), 0<f<°°, 

which has a pathwise unique, strong solution for t(-) given as Ti (y) := p l(o,e*>) (y) + 
a ] (y) by virtue of the NAKAO (1972) results, admits only a weak solution (and 
no strong solution) for T( ) equal to 

T 2 (y) = p l^oc) 0)- a l(_oo,o]()0- 

In both these cases the variance structure is the same, namely %\{y) = t\(y) — 

P 2 1(ohM + <7 2 1(-»,o]M- 

A rather similar situation obtains in Theorem 4.2 (and with obvious minor chan- 
ges, in Theorem 4.1 as well). To wit, both systems of equations (2.12), (2.13) and 
(2.14), (2.15) have the same covariance structure, namely s^{x\^xj) in (1.7), as the 
original system (1.1), (1.2); whereas the diffusion matrices for (2.12), (2.13) and 
(2.14), (2.15) are given respectively by 

e w (x UX2 ) = ( p ^ *W ) , l v (x u x 2 ) = (p^y a a \^A . 

\-P 1 {x 1 <x 2 } -^Hx 1 >x z }J \P 1 {x l <x 2 } GHx^}) 

(4.14) 

These two matrices differ only by changes of signs in the second row; yet the sys- 
tem (2.12), (2.13) of stochastic differential equations induced by the first of them is 
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strongly solvable, whereas the system (2.14), (2.15) induced by the second diffusion 
matrix is solvable only weakly. To paraphrase VERETENNIKOV (1982), page 448, in 
this latter case the diffusion generated by the system (2.14), (2.15) "cannot leave the 
diagonal {x\ —X2} strongly". 

We shall see in the next section that the original system (1.1), (1 .2), with diffusion 
matrix of the form 

L B (x h x 2 )=( pl ^ + a1 ^ Ql ° V (4.15) 

V P 1 {x l <x 2 }+ al {xi>x 2 } ) 

is also strongly solvable. A thorough discussion of weak and strong solutions, for 
all possible real square roots of the covariance matrix srf (x\,x%) in (1.7), appears in 
subsection 5.1. 

5 Strength, Restored 

In conjunction with the identity $( Xl > Xl \t) =3'C'iW(j) from (4. 1 3), the expressions 
in (3.12), (3.13) imply 

& B ^\t) C $ {X[ > X2) (t), V 0<f<°° (5.1) 

in the non-degenerate case p a > . It is an interesting question to settle, whether the 
reverse inclusion might also hold in (5.1), thus implying the strength of the solution 
to the system (1.1), (1.2) constructed in section 3. 

We do have such strength in the isotropic case p 2 = a 2 = 1/2 of equal vari- 
ances, as has been well known since VERETENNIKOV (1979, 1980, 1982). To obtain 
this property from first principles in our context, let us observe from (2.4)-(2.6) the 
representation W(t) = (B\(t) — 82(f))/v2 m me isotropic case, thus also the inclu- 
sion 

ft^fWcf'^W (5.2) 
for each t > 0; back into (2.5), (2.6) and in conjunction with (4.13), this gives 
£(x,,x 2 )( f ) = $(Wi,W 2 )^ c ^( B i*)(f), showing that (5.1) holds with equality in the 
isotropic case. 

More generally, the inclusion in (5.2) will follow whenever one is able to show 
that the modified TANAKA equation 

Y{t) =y-^J Q sgn(Y(s))ds+J (pl{y(,)>o}+tfl{rM<o})dBi(*) 

-J (pl{yW<Q} + t7l {J'W>o}) dB 2(*) ! 0<f<°o, (5.3) 

a consequence of (2.3)-(2.6), has a strong (that is, f?( B ^< B 2>— adapted) solution; then 
this strong solvability implies the inclusion 

jM^jISiA)^ v 0<f<°o, (5.4) 

that is, the strong solvability of the system (1.1), (1.2). Indeed, as Johannes Ruf sug- 
gests, the strong solvability of (5.3) implies the inclusion g^ 1 ^ (t) C $( Bl > B2 \t), 
and (5.4) follows then directly from the equations (1.1), (1.2). 
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Theorem 5.1. The system of stochastic differential equations (1.1), (1-2) admits a 
pathwise unique, strong solution; in particular, the inclusion (5.4) holds. 

Proof. We have seen already that the weak solution of (1.1), (1.2) constructed in 
section 3 is actually strong in the isotropic case p 2 = a 2 = 1/2 of equal variances; 
so we focus on the case p ^ a . 

From the results of Yamada & Watanabe (e.g., p. 3 10 in Karatzas & Shreve 
(1991)) and their "obverse" counterparts due to ENGELBERT (1991) and Chernyi 
(2001) we know that, in the presence of uniqueness in distribution, strong existence 
and pathwise uniqueness are equivalent; and in this case every solution is strong. We 
have already seen the that the solution of (1.1), (1.2) is unique in distribution, so 
we need only show that there is a strong solution; and by the argument before the 
statement of the Theorem, it is enough to prove that the solution of (5.3) is pathwise 
unique. Furthermore, it suffices to argue such pathwise uniqueness on an arbitrary 
but fixed time-horizon [0,T] of finite length T G (0,°°). 

A further reduction is that we need to prove such pathwise uniqueness only as 
regards the "driftless " modified TANAKA equation 

Y{t) =y + J Q (p 1 {yW>o} + ^l{y( 1 )<o})dj3 1 (s) 

o (pl{FM<o} + ffl{F(.s)>o})dj3 2 (i) (5.5) 

for < t < T , with j3 j (•) and J3 2 () independent, standard Brownian motions under 
a suitable equivalent probability measure Q; because then a Cameron-Martin- 
GlRSANOV change of measure brings us back to (5.3), with 

B i (t)=ft l (t) + —^, B 2 (t)=p 2 (t) + y^, 0<t<T 

independent Brownian motions under the original probability measure P . Let us ob- 
serve also, that uniqueness in distribution holds for the equation (5.5): every solution 
is a standard Brownian motion starting at F(0) = y, under the auxiliary, equivalent 
probability measure Q . 

Introducing the independent, standard Q— Brownian motions 

we can write (5.5) in the equivalent form 

Y(t) = y+ jTsgn(r(») dP(s) £±± *(0 . 
Pathwise uniqueness for this equation follows now directly from Theorem 1.1. □ 
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5.1 The Algebra and Geometry of Strength 

The matrices of (4.14), (4.15) are square roots of the covariance matrix sd (x\,X2) in 
(1.7). The general real square root of this covariance matrix is of the form 

.1/2/ v r , * (p 0\ (\ 0\ /coscp -sin<p\ 
^ ( X1 ,X2) = E{x u x2) = ^ a ) [q £ ) [ sin(p CQS(p ) 1 {X1 > X2} 

(o0\(l 0\ /cos#-sin#\ 
+ {0 p){0d) {sm$ cost? ) Hn^} - ^.b) 

where £ G { — 1, 1 } , 5 G { — 1,1}, < <p, # <2n, or equivalently 

I e/ 1/2 (x 1 ,x 2 ) = £(x 1; jc 2 ) = Z+l^^j+Z-l^^i (5.7) 

in the notation of (1.8). We also write the drift term of the generator (0.1) in vector 
form 

G(X U X 2 ) = l{ Xl < X2 }+ ( ^] 1{.V!>.V 2 } 

Each choice in (5.6) (or equivalently, (5.7)) leads to a system of stochastic differential 
equations with rank-based characteristics 

dX(f) = £/ 1/2 (X(»)dU(f)+G(X(f))df , (5.8) 

driven by a two-dimensional Brownian motion U(-) = (Ui(-), E/2G))' ■ This system 
admits a weak solution, which is unique in the sense of the probability distribution. 

Theorem 5.2. With the notation of (1.8) and ei = (1,0)', e 2 = (0, 1)', the (unique 
in distribution) weak solution of the system of stochastic differential equations (5.8) 
fails to be strong, if and only if 

(ei-e 2 )'2L = _(e 1 -e 2 )'Z + . (5.9) 

Proof. We note that the equation satisfied by the difference Y(-) =Xi(-) — X 2 (-) is 

dT(f) = (Cl 2 e2y ({E++E-) +sgn(T(0)(i: + -X-))dU(0-A S ga(Y(t))dt 

= -X sgn(T(f))df + ((ei - e 2 )' E + l {Y ( t )>0} + (ei -e 2 )'2L l{r(r)<o}) dU(f) 
= sga(Y(t))dM(t) + dN(t), 

(5.10) 

where we have set 

M(t):=<*^(E + -Z-)V(t)-kt and N(t) := (L + +Z_)V(t) . 

Elementary calculation shows that (M,N)(-) = 0, since (ei — e 2 )'Z± are two unit 
vectors on the plane. 

We observe next that, because the indicators l{x 1 (i)>x 2 (t)} an< ^ '^{x l (t)<Xi(t)} m ' e 
functions of the difference Y(t) =Xi(t) — X 2 (t), the cr-algebra ^^ Xl - X ^(t) is con- 
tained in the c-algebra g , ( f/ i f/ 2 7) (j) , for every t e [0,°°), by construction as in (5.8). 
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When the equation (5.10) for Y(-) is solvable strongly with respect to the planar 
Brownian motion U(-) = {Ui(-), Uii-))' ', we have the filtration comparisons 

5 M («) = 5 r Wc? M (r), thus c =3 lUlM2 \t) : 

the system (5.8) is then strongly solvable. When the equation (5.10) for Y(-) admits 
no strong solution with respect to U(-) , the system (5.8) admits no strong solution. 
There are now three possibilities. 

1. When (d — e?)'^+ = (ei — , the second line of the equation (5.10) for Y(-) 
simplifies to (2.3) which, as we have already remarked in subsection 2.3, has a 
pathwise unique, strong solution; so the system (5.8) is then strongly solvable. 

2. If both (ei — e2)'(X+±£_) are non-zero vectors, we can apply the Cameron- 
Martin-Girsanov theorem as in the proof of Theorem 5.1. To wit, we can 
restrict ourselves to a finite time horizon [0, T] and define a suitable equivalent 
probability measure Q under which N(-) and M(-) are continuous, strongly or- 
thogonal martingales. In this case the quadratic variations (M)(-) and (N)(-) are 
proportional, so the domination condition of Theorem 1.1 also holds. Then Y(-) 
is adapted to the filtration of (M(-),/V(-)), and this proves that the solution to (5.8) 
is strong in this case as well. 

3. Finally, it is possible that (ei — e2 )'(-£++£-) is zero; this condition is formulated 
in the statement as (5.9). Then the perturbation N(-) vanishes in (5.10), which 
becomes then an ordinary drifted Tanaka equation 

dY(t) = sgn(y(?))dM(f), M{t) = (ej - e 2 )'X+U(0 - It (5.11) 

driven by the Brownian motion M( ) with negative drift. From the Tanaka for- 
mula, the theory of the SKOROKHOD reflection problem, and (4.4), we note 

M(t) = \Y(t)\-\y\-2L Y (t), thus (t) = ffl y l (t) C , 0<r<-. 

We consider also the independent, one-dimensional standard Brownian motions 

U+(-) := (ei-e2)'X + U(-), := v'U(-), 

where the unit vectors E' + (e\ — 62) and V are orthogonal, and note S rU (f) — 
$ u +(t)V$ 3f {t) = ^{1)^3^(1), 0<t<°°. 

The equation (5.11) has no strong solution; that is, S rF (0 C 5 U (0 cannot possibly 
hold for < t < 00 . For if it did, then the process 

W(-) :=K(.)-y+A / sgn(T(f))df, 
Jo 

which generates exactly the same filtration as the process F(-) , would be a stan- 
dard one-dimensional Brownian motion and independent of 3f{-) , yet adapted to 
the filtration generated by U(-) ; but this is impossible, since 3 M (t ) = j' 7 ' (?) S 
$ Y (t) = y w (f) X ^(t) for < t < 00 would lead then to 

Z v (t) = d u+ {t)v^{t) - 3 M (t)v3 r (t) c $ w ( t ) v^ jr (f) c 5 U (0- 



Planar Diffusions with Rank-Based Characteristics 



21 



To summarize: the solution to the system (5.8) fails to be strong if, and only if, (5.9) 
holds. □ 

Remark 5.1. Under the condition (5.9), the diffusion vectors s + := (ei — ^2)' E + 
(for the half-plane {xi > jq}) and s_ := (ei — £2)' (for the half-plane {x\ < X2}) 
in the equation (5.10) point in exactly opposite directions. As a result, it becomes 
impossible for the planar diffusion X(-) = (X\ (-),X2(-)) of (5.8) to "escape strongly 
from the diagonal {x\ = x%Y'. 

Let us also note that s± are unit vectors, so the condition (5.9) holds if and only 
if their scalar product is — 1 . Working out this relation using the representation given 
in (1.8), we obtain the condition (1.9). 

Remark 5.2. Since (o 2 e - p 2 8) 2 + p 2 G 2 (l + s8) 2 = 1, the condition of (1.9) is 
equivalent to 

7T + iJ = <p + y (5.12) 

(modulo 2n), which involves only the angles # and <f> of (5.6), as well as the angle 
\jf G ( — 7t, 7i] determined from cos y/ = pc(l + e8) , sin y/ = o 2 s - p 2 8 . 

Example 5.1. The system of (1.1), (1.2), with diffusion matrix Eg in (4.15), corre- 
sponds to e = 5 = 1 , <p = # = and thus to y that has to satisfy cos \]f = 2 p <7 > , 
siny/ = <7 2 — p 2 ; there is no way that (5.12) can hold, so the system (5.8) is strongly 
solvable for all choices of p , a as in (1.3). 

Similarly, the system of (2.14), (2.15) with diffusion matrix E\y in (4.14) corre- 
sponds to e = —l,5 = l, <p = 0, # = — tz/2, and thus y/ = —n/2 as well. Once 
again there is no way for (5.12) to hold, so the system (5.8) is always strongly solv- 
able. 

Whereas the system of (2.12), (2.13) with diffusion matrix Ly in (4.14) corre- 
sponds to e = 1, 5 = — 1, <p = 0,# = —n/2, so y/ = n/2 for all choices of p , a 
as in (1.3). In this case (5.12) always holds, and the system (5.8) is never strongly 
solvable. 

Example 5.2. The matrix srf(x\,X2) in (1.7) has a total of 64 square roots of the 
form Ei l{ Xl> x 2 } +^2l{A,<.r 2 } with 

Among these, 48 lead to strongly solvable systems in the isotropic (p = a = 1 /y2) 
or degenerate (p a = 0) cases, whereas 56 choices lead to strongly solvable systems 
in all other cases. 

6 Joint Distributions 

The representations (2.22), (2.23) involve the triple (Y + (t),Y~ (t ),L Y (f)) as well as 
the independent random variable Q(t), where Y(-) is the diffusion process in (2.3), 
L Y (■) the local time of this process at the origin as in (1.6), and Q(-) an indepen- 
dent, standard Brownian motion. Thus, in order to compute the joint distribution of 
(Xi(t),X 2 (t)) via (2.22), (2.23), we need first to find that of (Y + (t),Y- (t),L Y (t)) . 
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In order to do this, we consider the "reference probability measure" P* , under 
which the process Y(-) becomes standard Brownian motion. According to (2.3) and 
the Girsanov theorem (e.g., Karatzas & Shreve (1991), section 3.5), we have 
for every t e [0,°°) the Radon-Nikodym derivative 



= exp(-A f sgn(Y(s))dY(s)-^- f sgn 2 (r(s))dsl 
3^(0 I Jo 2 Jo J 



exp A (\y\-\Y(t)\+2L ¥ (t] 



t\, (6.1) 



dP 
dPl 



thanks to the Tanaka formulas once again. Thus, for any Borel subsets A, B of 
[0, oo) and © of R, and with E+ denoting expectation with respect to the reference 
probability measure, we have 

P(F ± (r)eA, Y*(t) = 0, 2L Y (t) 6 B , Q{t) € 0) = expU\y\ - ^| X 



E 



exp {A (2L y (f ) - Y (t )) } l {y ± {t)eA _ yT (()=0 _ 2L v (f)6B} 



exp • 



%/27Tf 



(6.2) 



Remark 6.1. We have the classical result P*(Leb{0 < t < T : Y(t) = 0} > 0) = 
for the Lebesgue measure of the Brownian zero-set, thus P(Leb{0 <t<T: Y(t) — 
0} > 0) = holds from (6.1) for every T G (0,oo) . We deduce P(Leb{0 < t < °° : 
X l (t)=X 2 (t)} = 0) = 1. 

Let us also recall that the transition probability density function F(Y(t) € d^ | F(0) = 
y) = p r (v,^)d<^ for the one-dimensional "bang-bang" diffusion process F(-) of (2.3) 
is given by 



1 



V2wi 
for ^ > , and by 
1 



exp 



(v- 



-Ar) 



2r 



+ A e - 2 ^ f e-Wf'Wdu 
JyH 



(6.3) 



/27rr 



exp < 2Ay - 



It 



for 4 < . In particular, with y = the function 



(6.4) 



1 



'2w 



exp 



is evenly symmetric about the origin. Similar formulas hold for y < ; for the details 
of these computations, see Karatzas & Shreve (1984). 
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6.1 The Isotropic Case with y = x\-x2>0 

This equal variance case p = a = 1 /y2 affords the most straightforward computa- 
tion: from the representations of (2.25) and the independence of Y(t) and Q(t), we 
obtain 

F(x l{t ) e e d&) = exp (_ (iLtl^O! j d ^ 2 

2\f2nt L 2f J 

(6.6) 

for (§i,§2) € K 2 in the notation of (2.1), (6.3) and (6.5). The resulting transition 
probability density function 

^^=tyM- K,+ ^r v ' )2 } 

is continuous and strictly positive on all of R 2 , and of class <r # > °° on M 2 \ {(^1,^2) : 



6.2 The Degenerate Case with y = xj - x 2 > 

Let us focus now on the degenerate case with a = 0, thus p = 1 , and with x\ > X2 
as in Figure 1. From formulae (6.5.9)-(6.5.1 1), page 440 in KARATZAS & Shreve 
(1991), we have the joint probability distribution computations 



\(Y(t) Eda, 2L Y (t)edb 
a\+b+y) 



(|fl|+ ' + >' ) ^{- {H+t+y)2 \^b, aeR, b>0 (6.7) 



as well as 



Y(t)£da, 2Z/(f)=0 



1 ^xpf-^l-expl^M^ «><>, (6.8) 



V27Tf V I 2f J I 2f 

which are based on the theory of the so-called "elastic Brownian motion". Substitut- 
ing into (6.2) with = lwe obtain from these expressions 

¥[Y+(t) Gda, Y-(t)=Q, 2L Y (t)edb] = 
P[r(()edfl, y+(f)=o, 2L Y (t)edb] = 

fu ^1 \ A \\ ( a + b + y) J (a + b + y) 2 \ . „ 
= exp|A(y + fr- fl )- T rj- ^= -exp| j dad* 

= g -ax,. («+*+?) exp(- (a+fo+y ~ A?)2 Udfc, a>0,*>0. (6.9) 
V27Tf 3 I 2? J 
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In a similar fashion, we obtain 

P(F ± (f)eda, Y*(t) = Q, 2L Y (t)=0) = 



2 J 72lr7 V I 2f J I 2f 

1 ^p/-^^)- e --exp/-to^)^ (6.10) 



\/27r7 V I 2f J ^ { 2t 

for a > . This expression vanishes, as it should, for y — : in this case the process 
Y(-) starts accumulating local time at the origin straightaway, that is, P(L Y (t) > 0) = 
1 holds for every t £ (0,°°) . 

• We set out to compute the joint distribution of the random vector (Xi (f ),X 2 (f )) 
for given, fixed t > 0. From Remark 6.1 and either (2.26) or (4.8), it is clear that this 
distribution is supported in the planar region *Bi U 5$2, the union of the two blunt 
(135° -) wedges 

Bi := { (6 , &) : l\ > & , & < *2 +gt} , B 2 := { (6, &) : & < & , & < Jd +f f } . 

• Let us consider the wedge 03 1 first. With given real numbers <jji , ^2 that satisfy 
£l >fe, £2 <X2+gt, setting 

£:= {(fl^)e(O,-) 2 ; X2 + g t+a-(b/2) xa+gt- (ft/2) <&}, 

and with the help of the expressions (2.26) and (6.9), we obtain 

p(Xi(0>§i,x 2 (0<&) = 

( a + fc +^) exp (J a + fc +/-^) 2 \ dad , = 



\/2^ I 2f 

r fr ex P (- (a+fe+ ^ Af)2 ldfl)d.. 

We differentiate this expression, first with respect to § 2 , then with respect to §1 ; 
recalling the notation of (2.1), we obtain 

P(Xi(f)ed§i,A2(t)ed&) = 

2e -n S - M . l.-3fe+y2 g , .„ p {_J_ (| ,_ 3;2 + z + v ^} d; , d|;; 

|l>&, fe<*2+gf- (6.11) 

On the other hand, with §1 > ^2 = ^2 + gf there is no accumulation of local time at 
the origin over the interval [0,f] , so the expressions of (2.26) and (6.10) give 

p(z,(0 e J*(0 = &) = ■ (ex P {- (Ilz^!) 2 

(6.12) 
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• Next, we consider the wedge 2$2 ; equivalently, we work on the event {X\ (t) < 
Xlit)} , on which the initial order x\ > X2 stands reversed at time t and L Y (f) > 
holds a.e. In particular, the joint distribution of the random vector (X\(t),X2(t)) 
assigns zero mass to the region {(^1,^2) : §2 > %\ = *2 + gt}, as we have already 
observed. 

With given real numbers §1 , £,2 that satisfy §1 < %\ < X2 + gt, denoting 
£ := {{a,b)e {0,oo) 2 :x2+gt-(b/2)<% u x 2 +gt +a - (b/2) > fe} , 
and with the help of the expressions (2.26) and (6.9), we obtain then 



»(Xi(0<Si,X2(*)>S 



2 



{a + b + y)_<Ja + b + y-kt? 



2{ X2 +gt-{,i) \JiL+&-x2-gt V2nt 3 I 2f 

Differentiating this expression, first with respect to <jji and then with respect to 
we obtain the following expression for the probability density function: 



2r"^,). fe- 3 ^i+^ .exp(-i-fe-3^+ Z + vr)-U^dfc; 



"(x l {t)ed^ u x 2 (t)ed^ 

v / 2^ I 2f 

§i<&, £i<x 2 +gf. (6.13) 

This is the same as the expression on the right-hand side of (6.1 1), except %\ , £2 
have now traded places. 



The following Figure plots the joint Probability Density Function of (Xi (t), Xj (f)) . 




Fig. 2 Joint density: g = h 



= 1,0" 



= 0, p = 1, jci =x 2 = 0, t = 1 (left), f = 2 (right). 
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Remark 6.2. The joint distribution for the ranks R\(t) =max(X\(-),X2(t)) , Riit) = 
mm(X\(t),X2(t)) is supported by the planar region {(p\,pi) ■ Pi > p2, P2 < ^2 + 
gt}; from (6. 1 1)-(6. 13), or directly from R 2 (t) = r 2 +gt-L Y (t) and Ri(t) = r 2 + 
gt + \Y(t)\ -L Y (t) and (6.7)-(6.8), it is computed as 

F(R l (t)ed Pl ,R 2 (t)edp 2 ) = 

4(2gt+z + Pi-3po) f w s 1 / \ 2 1 

— * ^Ynti ^ •exp|-2^(p 1 -p 2 )- — (^z + vf+pi-3p 2 J |dpidp 2 ; 

Pi <x 2 +gt, pi > p 2 (6.14) 

and 

P(*i(0edp 1 ,* 2 (0 = P2) = 

Pi >P2 =x 2 +gt . (6.15) 

Remark 6.3. Once again, the probabilities in (6.12) and (6.15) vanish for y = 0, 
i.e., when the two particles start at the same point. In this case the distribution of 
(X\(t),X2(t)) is absolutely continuous with respect to Lebesgue measure on the 
plane, with probability density function (<jji , §2) given by (6. 1 1) on E,\ > £ 2 , £2 < 
X2 +gt (in the wedge 2*i), and by (6.13) on %\ < £ 2 , 4i < x 2 +gt (in the wedge 
B2). 

Even in this case, though, there is a discontinuity along the front {(§1 , £2) S : 
%i^t;2=x + gt} (cf . Figure 2), of size 

S L -{-"l«'-fcl-F0«'-«'l- A ') 1 }- 



6.3 The Non-Degenerate Case with Unequal Variances 

More generally, and by virtue of (4.2), (6.9) and (6.10), the quadrivariate joint distri- 
bution of (F+(?),y-(f),L F (f), Q{t)) is given as 

P[T ± (f)edfl, Y T (t)=Q, 2L Y {t)edb, Q(t)£d&] = fi(a,b,-&)dadbdi} (6.16) 

fora>0, b>0, #eM,and 

P[F ± (f)eda, 7*0) =0, 2L y (»=0, Q(» € d#] = / 2 (a, 0) dadtf (6.17) 

fora>0, #eM,with 

/i(a,6,dj=e ■ 27rf2 — exp I — j, (6.18) 



Planar Diffusions with Rank-Based Characteristics 



27 



h{a,-d)= — exp f-exp<^ -2Aa 



2nt V I 2f J I 2f ) J 

(6.19) 

In conjunction with the skew representations of (2.22) and (2.23), now written in the 
form 

Xi(t) = xi+nt + %(t), X 2 (t) =x 2 + llt + %,{t) 

with 

%(t) := p 2 (F+(f) -y+) - a 2 (Y-(t)-y-) - yL Y (t) + poQ{t) , (6.20) 
f 2 (r) := -a 2 (y + (f)-3' + )+p 2 (^(f)-3' _ )-7i F (0+P^e(0, (6-21) 

it clearly suffices to compute the joint distribution of (W\ (t ), . This is facilitated 

by the observation that the system of (6.20), (6.21) can be "inverted", in the sense 



Y + {t)~y + = - (/ % (t) + & W)) + L * (0 - Y fi(0 • 



Y-(t)-y- = - [a* K<t)+p*W))+L'(t)-!j- Q{t) . 

To proceed further, the cases y > , y < and 7 > , 7 < have to be considered 
separately. We shall discuss briefly only the case y = x\ — x% > 0, 7 > (i.e., p > 
a > 0) In this case, the joint probability density function 

P(f 1 (f)edv/i, f 2 (f)edv/ 2 ) = Vf(Vi.Vfe)dyidv<>2 

of (<Fi(/), «f*>(r)) in (6-20), (6.21) is given by 

7<P,(y/i, 1/2) = 

•0(Vl,V2) 



/i(6-b*(vri,%,*))d6+/a(-b*(Vi,Vi.*))|di> 

/ 1 (/7-b*(v/i,r2,^))d^dt>, (6.22) 



where we have set 

:= yi+ 2 pa +ry ' & *^,V2,0) : = ^(*"*(^.^)) ■ 
After some calculations, (6.22) reduces to 

% (Yi, ¥2) = -j= (1 + * (A+ ( Vi , ^2)) exp I - i ( v/i + V2 + At? ) 2 J 

- e 2 ^'0(A_(v/i ! r2))exp{-i(v/i+V/2+27y + A7f) 2 }] , (6.23) 

where <£>(■) := (1/s/2k ) e ~*V 2 dz is the cumulative standard normal distribu- 
tion function and 

y 

vTA ± (y/i,y/2) := — (Vi +V / 2 + 7y)±2pcr(y=FAf) , (ift,ifc)eK ■ 
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7 Time Reversal 

Let us consider now the time-reversed versions ("free" and "anchored", respectively) 

Y(t) := Y(T-t), W(t) := W(T-t)-W(T), 0<t<T, (7.1) 

over a given time horizon [0, T] of finite length, of the diffusion process Y(-) in (2.3) 
and of the Brownian motion W(-) which drives that equation. Both of the processes 
introduced in (7.1) are adapted to the backwards filtration F = {J^(f)}o<?<r defined 
as 

&(t):=o(Y(T))V3*(t), $ W (t):=o(W(6);0<8<t). (7.2) 

Note that the process W(-) is Brownian motion with respect to the filtration F w = 
{3 W (t)}o<t<T ; this is because W(-) has continuous paths, and is easily checked to be 
an F w — martingale with the same quadratic variation as Brownian motion. However, 
W(-) is only a semimartingale with respect to the larger backwards filtration F; its 
semimartingale decomposition is provided by the fact that the process 

W # (t) := W(t)- J' q(T -s,Y{s))ds = W{T -t) -W(T) - J q(s,Y(s))ds 

(7.3) 

for < t < T , is a standard F— Brownian motion. Here 

q(T,|) :=~logp T (y,S), (t,§)6(0,«)xR (7.4) 

is the logarithmic derivative of the transition probability density function in (6.3), 
(6.4). 

Remark 7.1. This result is proved as in Meyer (1994), who notes also the following 
corollary: the (once-more-time-reversed) process 

J7(f) :- W # (T -t) - W # (T) - Y{t)+f (q(*,^(*)) +^ sgn(F(j))) ds 

= W(t) + J\(s,Y(s))ds = J' (sgn(F(>))dV b (s) + q(s,FOO)ds) , 0<t<T 

(7.5) 

is also Brownian motion with respect to its own filtration, and is adapted to the filtra- 
tion F F = F w (though not to the filtration F v ; recall the filtration identity (4.4) in 
this regard). 

The process »)(•) of (7.5) is independent of the random variable Y(T). In fact, for 
every given t E [0,T), the a algebra #1 (t ) = <j(t] (s) , < s < t) is independent of 
the (J— algebra o(Y(9) , t < 9 < T) ; in particular, of the random variable Y(t). This 
shows that the inclusion ^(f) C "S Y (t), quite obvious from (7.5), is strict; or, put 
another way, that the stochastic integral equation 

Y(t)=y-J t (Xsga(Y( S ))+q(s,Y(s))^ds + ri(t), 0<t<T (7.6) 

cannot possibly have a strong solution. □ 
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A result by HAUSSMANN & Pardoux (1986) (see also the expository papers 
by Pardoux (1986) and Meyer (1994), among many other works on this subject) 
is that, with respect to the backwards filtration F = {,^(t )}o<r<r , the time-reversed 
process Y(-) in (7.1) is not just a semimartingale but also a diffusion process driven 
by the F— Brownian motion in (7.3), of the form 

Y{t) = ?(0)+ [b (T -s,Y(s)) ds + W # (f), 0<r<7\ (7.7) 



Here, the new (backward) drift function b( - , •) is given in the notation of (7.4) by the 
generalized NELSON equation 

b(T,|) = Asgn(4)+q(T,4), (T,§)G(0,r]xR. (7.8) 
Remark 7.2. In the special case y = 0, the function of (7.4) takes the form 

■ ^(2X-^(p^(t,-^)+2X 2 e 2 ^ J^(p^(t-u)du^j , 4<0, (7.9) 
and q(T,|) = -q(T,-|) for % > 0, T > with 



<pW(T^):=-^e X p{-ii±^^. (7.10) 
In particular, we obtain from the equations (7.9)-(7.10) the explicit expression 



t>(T,£) = A S gn(|) - \cpW(z, |4 1) +A e - 2A ^ J <pW(T,-u)du^ 

'(2A + ^)<pW(T,|4l)+2A 2 e - 2 ^l |" -«)d«) , 4 ' 



(7.11) 

The singularity at z = of the drift in (7.1 1) is of the "bridge" type: it ensures that, 
as the time-to-go T I decreases to zero, the backward diffusion Y(-) zooms into 
the prescribed terminal condition Y(T) — (the initial condition Y(Q) = of the 
forward process), as it must. □ 

Let us also note that the semimartigale local time L Y (•) of the backward diffusion 
Y(") , and the semimartigale local time L Y (•) of the forward diffusion Y(-) , are linked 
via 

L Y {t) = L Y (T)-L Y (T -t), 0<t<T. (7.12) 



30 



E. Robert Fernholz et al. 



7.1 A Time Reversal for Names 

We consider now the "free" time-reversals 

Xi(t) := Xi(T-t), X 2 (t) := X 2 (T-t) , 0<t<T (7.13) 

of the components of the vector process (Xi(-),X2(-)) constructed in section 3, as 
well as their "anchored" versions 

Xj(t):=Xj(t)-Xj(0)=Xj(T-t)-Xj(T); < t < T, j = 1, 2. (7.14) 

By analogy with (7.2), we also look at the new backwards filtration F = {g(f )}o<r<r 
given by 

ff(0:=ff(r(r))v$<&^(0, with ^■ iv) (t)^a(Q(e),w(e);o<e<t) 

generated by the random variable Y(T) and by the time-reversed versions 

Q(t) = Q(T-t)-Q(T), W(t)=W(T-t)-W{T), 0<t<T (7.15) 

of the independent Brownian motions W(-), Q(-) of (2.19), (2.21) in the manner of 
(7.1). In particular, Q(-) is independent of both W(-) and Y(T), thus also of the 
Brownian motion W # (-) in (7.3). 

Then the skew representations (2.22) and (2.23), along with the notation of (7.1) 
and the local time identity (7.12), imply that the "anchored time-reversals" of (7.14) 
are F— adapted and are given by 

Xi (t ) = -\it + p 2 (r+(r) - ?+(0)) - a 2 (?- (t) - Y- (0)) + yL ? (t) + pa Q(t) 
= -Ht + J (p 2 l { ? W >o} + ° 2 1 {?(s)<0}) AY{s)+2yL Y {t) + poQ{t) , (7.16) 

and 

X 2 (t) = -fit-a 2 (Y+(t) -?+(0)) +p 2 (Y-(t) - Y-(0))+ Y L 9 (t) + poQ{t) 

= -Lit-J^(p 2 1 {?(J) < 0} + <? 2 1 { ? (J)>0} ) dY(s)+2rL ? (t) + poQ{t) , (7.17) 

respectively, thanks once again to the Tanaka formulas. 

We recall now the "backwards dynamics" of (7.7) as well as the notation of (7.4), 
(7.8), and write these equations in the time-reversed skew representation form 

X l {t)=X l {t)-X l {0) = 

= f ( h hx,( s )>x 2 ( s )}-S\2 l{s) <x 2 ( s) ^ ds + (P 2 - o 2 )l}^\t) + poQ(t) + 

j£(p 2l { f lW >f 2W} +^^ 

(7.18) 
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and 

X 2 (t)=X 2 {t)-X 2 {0) = 

(7.19) 

Remark 7.3. A somewhat interesting dichotomy emerges. In the case of equal vari- 
ances (p 2 = a 2 = 1/2) these anchored time-reversals are F— adapted Ito processes: 
the bounded variation terms in their semimartingale decompositions are absolutely 
continuous with respect to LEBESGUE measure. In the case of unequal variances 
p 2 7^ (7 2 , terms which are singular with respect to LEBESGUE measure appear, and 
are governed by local time. 

To the best of our knowledge, this is the first time such a feature is observed in the 
context of time-reversal of a "purely forward" stochastic differential equation with- 
out reflection; its occurrence and significance need to be understood further. For a 
similar but different phenomenon, in the context of time reversal of Brownian mo- 
tion reflected on an independent time-reversed Brownian motion, see SOUCALIUC 
& Werner (2002) (as well as Soucaliuc et al (2000), Burdzy & Nualart 
(2002)). 

7.2 A Time Reversal for Ranks 

By analogy with (7.13), we introduce the time-reversed versions 

*i(f) := Ri{T-t) = max(Xi(0,X 2 (0)> R 2 (t) := R 2 (T-t) = min(Xi(r),X 2 (0) 

(7.20) 

of the ranked processes in (4.1) for < t < T . We have R\ (t) +R 2 (t) =X\ (t) +X 2 (t) 
and Ri(t) —R 2 (t) — \X\(t) — X 2 (t)\ — \Y(t)\, so the time-reversed skew representa- 
tions of (7.16), (7.17) cast the "anchored" versions of the processes of (7.20) in the 
form 

Ri(t) :=/?i(0-£i(0) = -lit + p 2 (\Y(t)\~\Y(0)\) + Y L 9 (t) + p aQ(t) , (7.21) 

R 2 (t) :=R 2 (t)-R 2 (0) = -iit-o 2 (\Y{t)\-\Y{Q)\)+ylJ(t)+poQ{t). (7.22) 

In conjunction with the reverse-time dynamics of (7.7), (7.8), the Tanaka-Meyer 
formulas give now 

|?(0|-|?(0)| = Af + J^sgn(Y{s))q(T-s,Y{s))ds + V # (t)+2L f {t), (7.23) 
where 

V # {t) := / sgn(F(s))dW # (», 0<t<T (7.24) 
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is a Brownian motion of the backwards filtration F , and is independent of the Brow- 
nian motion Q(-) . 

Substituting the expression of (7.23) back into (7.21) and (7.22), and recalling 
(2.1), (2.24), we obtain the dynamics 

£i(r)-£i(0) = ht + p 2 J' sgn(Y{ s ))q(T-s,Y( s ))ds + P V?{t) + (4p 2 -l)L/(t), 

(7.25) 

R 2 (t)-R 2 (0) = -gt-o 2 f sgn (Y(s)) q(T -s,Y(s)) ds + aV 2 # (t)~ (4a 2 -l) L 9 (t) . 

(7.26) 

Here the processes 

Vf(.) := pV # (-) + oQ(-), V 2 %) := p Q(-) - oV # (-) (7.27) 

are independent, standard Brownian motions of the backwards filtration F . 

Comparing the equations of (7.25), (7.26) with those of (4.5), (4.6), we see that 
V\ (•), V* (•) play in the context of time reversal the same roles that the Brownian 
motions Vi(-), V2O) play on the forward context: to wit, that of independent, stan- 
dard Brownian motions associated with individual ranks. 

7.3 Steady State 

Finally, let us note that the diffusion process Y(-) of (2.3) has invariant distribution 
with double exponential probability density function 

p(§) = A«- n l*l, i;eR. 

For this function, the analogue q(§) = (d/d£A logp(§) of the logarithmic derivative 
in (7.4) becomes q(^) = — 2A sgn(^) in the notation of (6.1), and from the general- 
ized Nelson equation (7.8) the drift in the backward equation (7.7) becomes 

m = -X sgn(|). 

This reflects the fact that the diffusion process Y(-) of (2.3) is strictly time-reversible 
when started at its invariant distribution, that is, the processes Y(-) and Y(-) are 
then identically distributed; and then the Brownian motion in (7.5) takes the form 
tj (t ) = J Q ? sgn (Y(s)) (dV b (s) - 2X ds) , < t < T . 

In such a setting, the equations of (7.6) and (7.18), (7.19) continue to hold, now 
with q(-) = — 2X sgn(-); whereas the equations of (7.25), (7.26) assume the rather 
concrete form 

Ri(t)-*i(0) = (h-2(g + h)p 2 )t + P V*(t) + L S ^(t), 

R 2 (t)-R 2 (Q) = (2 (g+h) a 2 -g)t + aVl{t) - L«^(t) , 

where L R ^ Rl ( ) = 2L Y (-) is the local time accumulated at the origin by R-i(-) — 
R 2 (-) = \Y(-)\. 
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8 A Generalization of the perturbed Tanaka equation 

We have the following generalization of Theorem 1.1, which partially answers a 
question posed by Professor Marc Yor. In what follows, we shall agree to denote 
by (£?)(■) the quadratic variation of the continuous local martingale Y(-) 

in the decomposition of a continuous semimartingale 3?( ) = 2 (0) + "f (■)+£§(■) , 
where 3§( ) is a continuous, adapted process of finite variation on compact intervals. 

Theorem 8.1. Let f : R — > R be a function of finite variation, that is, f = f + — /_ 
where f± : R — > R are increasing. On some filtered probability space (Q, J?",P), 
F = {$(t)} t >Q, consider two continuous local martingales M(-), N(-) which satisfy 
the conditions of Theorem 1.1, and a continuous, adapted process A(-) with A(0) = 
and finite total variation A(f) on compact intervals of the form [0,?]. 

Then pathwise uniqueness holds for the stochastic differential equation 

Y(t) =y + J f(Y(s))dM(s)+A(t)+N(t), 0<t<™. (8.1) 

Proof. We shall prove the statement for / = /+ — /- with bounded, increasing f±. 
The general case will then follow by considering /'"' = —/I"', where f± ~ 
min (n,max(— n,/±)) , n G N are truncated versions of /. If pathwise uniqueness 
holds in the bounded case, then the solution is pathwise unique up to the stopping time 
sup#„, where = inf{f > : |/±(T(f))| > n}. That is, if we have two solutions, 
X(-), Y(-) to (8.1) on the same probability space, issued from the same initial value 
and driven by the processes M(-), Af(-) and A(-), then X(t) — Y(t) holds for t < 
i^(X) A £(F), where £(Y) denote the lifetimes of the processes X(-) and T(-), 

respectively. 

Suppose then that the processes X(-) , Y (•) both solve the equation (8.1), so their 
difference 

D(t) :=X(t)-Y(t) = £(f(X(s))-f(Y(s)))dM(s), < t < ~ (8.2) 
is a continuous local martingale with quadratic variation 

(D)(t)= Jjf(X(s))-f(Y(s))) 2 d(M)(s)< \\f\\ TV Jjf(X(s))-f(Y(s))\d(M)(s) 

(8.3) 

Here jj/jlj-y is the total variation of / over the real line, and we have used the ele- 
mentary comparison (f(x)-f(y)) 2 < ll/Hrv \f{x)-f(y)\, V (x,y) e R 2 . Following 
Le Gall (1983), we shall establish the estimate 

E Io ^dUt ~ c ' W f Wrv ' ™V E ( 2L(M) ( r ' fl )) (84) 

«6[0,1] 

where, for each u g [0, 1] , the quantity lS u \T,a) = A z< "' (T,a) is the local time ac- 
cumulated at the site a E R during the time interval [0, T] by the continuous semi- 
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martingale 

Z«(.) := {\-u)X(-) + uY{-) = 

y + J o ((l-u)f(X(t)) + uf(Y(t)))dM(t)+A(-)+N(-). (8.5) 

To this end, we introduce a sequence {fk}keN C ^ (M) of continuous and con- 
tinuously differentiable functions that converge to / pointwise, and are bounded in 
total variation norm, that is sup£ GN ||A||rv < °°- Since limsup^ ||/t|| r y < ||/||j-y orJ " 
viously holds, it is only possible if / is of bounded variation and in this case an 
approximating sequence is easily obtained, e.g., by mollifiers. We note the identity 
fk(X(t))-f k (Y(t)) = (X(t)-Y{t)) /J f k (zW(t)) du, as well as the comparison 



\f k (X{t))-f k (Y{t))\ 

x ^_Y( t ) 1 {x(t)-Y(t)>S}<i(M}{t) 

E [ T \f k (Z (H) (0)| l{x W -r W >5} d(M) (0 ) du 



E 



for 5 > 0. We also note that the orthogonality of M(-) and N{-) implies (zW)(-) > 
(N)(-) in conjunction with (8.5). Let us recall now the domination condition (1.11), 
which gives 

(M)(-)<c(N)(-)<c(zW)(.)<C(N)(-) 
for a suitable real constant C > c. We deduce from this 

-1 / i-T 



jf [eJ o \fi(zW(t))\l {D{t>S] d{M)(t))du 
<c f Uj^\fl(Z^{t))\d{Z^){t)\du 



E / \fj c (a)\2L (u] {T,a)da)du 



<c- sup e(2L w (7»V / da < c-||/*|| rv - sup E (2L {,l) {T,a)) . 



I TV 

am V 7 JR aeE 

u€[o,l] "£[0,1] 



Letting k j" °° and then 5 1 , these estimates give 



E 



J§T^0 ^('KO^^W < c -ll/llrV sup E(2&\T,a)) 

»e[o,l] 



and, in conjunction with (8.3), the claim (8.4) as well. 
Suppose now we can show 

r(«). 



a EH 
»£[0,1] 



sup E[2L w (T,a)) < °°. (8.6) 
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On the strength of (8.4), this then implies E/ r (D(f l[o(t)>0} d(D) (f ) < °° for 
each T € (0,°°) ; and arguing as in Le Gall (1983), Lemma 1.0 (see also Exercise 
3.7.12, pages 225-226 in Karatzas & Shreve (1991)), we deduce that the local 
time L(-) = A D (- ,0), accumulated at the origin by the continuous local martingale 
D(-) of (8.2), is identically equal to zero. But then, by the Tanaka formula once 
again, we obtain that \D(-) \ is a local martingale, thus also a (nonnegative) continuous 
supermartingale with |D(0)| =0, and consequently X(-) — Y(-) =D(-) = 0; that is, 
pathwise uniqueness holds. 

• The property (8.6) is checked by standard methods: with the help of the Tanaka 
formula 

|ZM (T)-a\ = |ZW (0) - a \ + J sgn (Z^ (?) - a) dZ^ (f ) + 2 L<"> (T, a) 
and the Ito isometry (e.g., KARATZAS & Shreve (1991), p. 144) we get 

e(2zM(7») <e zW(r)-zW(o) +E 1 / 2 ((z<"))(r)) +E(A(r)) 
<2E 1 / 2 ((zM)(r)) +2E(A(r)) 

< 2 VC E 1 / 2 ((N) (T)) +2E (A(T)) . 

This last quantity does not depend on a € R or u £ [0, 1] ; if it is also finite, we are 
done. 

If not, we deploy standard localization arguments: to wit, we consider the stop- 
ping times 

T,„ = inf { t > : max (A(t),(N) (?)) > m} , 

and deduce D( • A T m ) = from the above analysis, for every m € N, P— a.s. But 
then lim m ^«, T,„ = °° also holds P-a.s., and this leads to D( ■ ) = once again. □ 

Our next result covers cases discussed in Theorems 5.1 and 5.2. More importantly, 
it generalizes Theorem 8.1, by replacing in the equation (1.10) (respectively, in the 
equation (8.1)) both driving local martingales M(-) and N(-) by semimartingales. 
This generalization can be construed as an analogue of the results in ZVONKIN (1974) 
and VERETENNIKOV (1979). With bounded, measurable / : W -> R m , these authors 
show (for m = 1 and for general m € N, respectively) that, even in situations in 
which the ordinary differential equation d¥(t) = f(Y(t))dt might not be solvable, 
the addition of a Brownian perturbation W(-) as in 

dY(t) = f(Y(t))dt + dW(t) 

restores to the differential equation a pathwise unique, strong solution. 

Proposition 8.1. With the same assumptions and notation as in Theorem 8.1, and 
with r(-) a continuous, adapted process of finite first variation on compact intervals, 
pathwise uniqueness holds for the stochastic differential equation 

Y(t) =y + J' f(Y(s))d(M( S )+r( S ))+A(t)+N(t), 0<r<~ (8.7) 

provided that either 
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(i) the process !"(•) is increasing, and the function /(•) is decreasing; or 

(ii) the process F(-) is decreasing, and the function /(•) is increasing; or 
(Hi) the process r( ) is of the form 

r(-)= [ <p(t)d{M)(t) 

Jo 

for some progressively measurable process (p(-) which satisfies the integrability 
condition f„ \(p(t)\ 2 d(M)(t) < oo, V re(0,°°). 

Proof. We show exactly as before that the continuous semimartingale D(-) — X(-) — 
y(-) accumulates zero local time L(-) = A D ( - ,0) = at the origin, so we have 

\D(-)\=J o n(t)(dM(t)+dT(t)), (8.8) 

with 

Tl(t):=sgn(X(!)-Y(t))(f(X(t))-f(Y(!)))- 

Under either of the conditions (7) or (7/), the process |D(-)| is now a continuous 
local supermartingale thanks to the representation (8.8), thus a true supermartingale 
(by Fatou's lemma) since it is nonnegative. As we have D(0) — 0, we conclude 
D(-) = just as before. 

Under the condition (Hi), we note first that the process <p( ) satisfies also the 
integrability condition f Q T \(p(t)\d(M)(t) < °o for every T e (0,°°) , by the KUNITA- 
Watanabe inequality; cf. Karatzas & Shreve (1991), p. 142. We introduce the 
continuous local martingale 

K(-):= [ (p(t)dM(t), note r(-) = (M,K)(-) , (8.9) 
Jo 

and summon the "stochastic exponential" 

30) :=« P j-*(.)-i <*>(•)}, 

the unique solution of the "simplest stochastic integral equation" (in the terminology 
ofMcKEAN(1969)) 3(-) = 1 -J 3(t)dK(t). 

We want to show D(-) = 0. Observe that IT(-) := 3(0 \D\ (•) is a local martingale 
by the product rule of the stochastic calculus: 

n(-) = [j(t)d\D(t)\+ f\D(t)\d3(t)- [n(t)3(t)W,m) 

Jo Jo Jo 

= /' 3(0*7 «dM( f )+ /'|D(r)|d3(0- 
Jo Jo 

Here M( ) and 3(-) are continuous local martingales and the integrands are locally 
bounded, hence both terms on the right are continuous local martingales. Now, JT(-) 
is a nonnegative, continuous local martingale, starting from the origin at time zero, 
so it must stay at the origin at all times; and since 3(0 is strictly positive, this leads 
to the desired conclusion D(-) = 0. □ 
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8.1 Maximality 

Let us place ourselves in the non-degenerate case p a > and in the constructive 
setup (synthesis) of section 3, where we work with the filtration J?( Wl > Wl ' generated 
by the planar Brownian motion (Wi (•) , W 2 (•) ) . 

In the terminology of BROSSARD & LEURIDAN (2008), we shall say that a one- 
dimensional Brownian motion /?(•) in the filtration F^ 1 '"'^ , is 

(i) complenentable in F ( Wl , if there exists an independent one-dimensional Brow- 
nian motion Tj(-) in the filtration F^' W ^ , and for which F^'^ = ; 

(ii) maximal in F^ Wl ' Wl \ if for any one-dimensional Brownian motion t?( ) in the 
filtration F^'^ and for which $P(t) C $*(t) holds for all < t < °o, we 
have $°(t)=3 P (t), 0<f <°c. 

BROSSARD & LEURIDAN (2008) show that complementability implies maximality; 
it is still an open question whether the reverse is true. 

In section 3 we start with the planar Brownian motion (W\ (■), W 2 (-)) and con- 
struct 12 one-dimensional Brownian motions in the filtration F' Wl ,W2 ' . We shall ex- 
amine these properties for each of them. 

Let us start with the independent one-dimensional Brownian motions B\(-) of 
(3.12) and B 2 (-) of (3.13). In view of (5.1), (5.4), Theorem 5.1 and Proposition 4.2, 
we have then ^ B ^\t)=^ x ^\t)=^ w ^{t), 0<t <°o.Thus, each of Bi(-) 
of (3.12) and B 2 (-) is complementable (by the other one), therefore also maximal, in 
the filtration F^ w '' W i\ 

It follows also fairly directly from (3.1) and (3.2), that the independent Brownian 
motions W(-) and U (•) complement each other in F' Wl,W2 ); the same is true of 
U (•) and W (•) . Thus, all these one-dimensional Brownian motions are maximal in 

On the other hand, let us consider the one-dimensional Brownian motion V (•) of 
(3.5). This process is adapted to the filtration F w generated by the Brownian motion 
W(-) of (3.1); this is aconsequence of the representation V\-) = J Q sgn (Y(t)) dW(t) 
in (2.17), and of the strong solvability of the stochastic equation (1.4). But we also 
have $ v \t) = £l F l(?) C ^(f) = $ w (t), < r < 00 from (4.4), so V b (-) cannot 
possibly be maximal in the two-dimensional filtration F^ 1 ™ . 

As a corollary of these considerations, we observe also that a linear combination 
of independent Brownian motions which are maximal, such as V\-) = p Vi(-) — 
crV2(-) in (3.5), can fail to be maximal. 

As for the Brownian motion Q "(•) of (3.4), it can fail to be maximal. Indeed, 
it is clear from (3.4), (3.5) and the observations in the previous paragraph that in 
the isotropic case p = <J this Brownian motion Q (•) = V (•) is not maximal in 
F (Wi,w 2 ) It turns out that q\>^ is max i ma i i n pW*) , however, if p ^ a. To see 

this, recall the notation 8 = 2pa and 7= p 2 — a 2 from (2.18), and observe 

8 f sgn(Y{s))dQ\s) + YW\t) = {8a + yp)W l (t) + {Sp -ya)W 2 (t) = 
Jo 

= pWi(t) + aW 2 (t) = W(t) = Y{t)~y-X [ sgn(r(s))ds; 

Jo 
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equivalently, 

Y(t) = y + J^sgn{Y(s))d(8Q b {s)+?.s)+yW\t), 0<t <~. 

Then, since we are assuming 8 = 2pa > and y = p 2 — a 2 ^ 0, it follows from 
Proposition 8. 1 (iii) that pathwise uniqueness holds for the above equation; and since 
the equation admits a weak solution, this solution is actually strong. In particular, 
the process F(-) is adapted to JyQ ,w \ and we arrive at the filtration identities 
F (e b .w^) = F (Qt,wKY) = F (Qt,wKw) = F (W U W 2 ) Xhus com pie m ents the Brown- 

ian motion Q b {-) in F< w i^). 

Recall now the one-dimensional Brownian motion Q(-) of (3.5), which is inde- 
pendent of V b (-); it is also independent of the one-dimensional Brownian motion 
W(-), and indeed we have $Wi> w d(t) =$( w <®(t), < t < <*> from Proposition 4.2. 
In other words, the one-dimensional Brownian motion Q(-) is complementable (by 
W(-)), therefore also maximal, in the two-dimensional filtration F( Wl ' Vl ' 2 ' . 

Let us consider now the one-dimensional Brownian motions Vi(-) and V2(') 
of (3.3); we claim that they are complementable (by W2(-) and W\(-), respectively; 
though not by one another!), therefore also maximal, in the two-dimensional filtration 
p(Wi,w 2 ) jo see the first claim (the second one is argued in a completely analogous 
manner), let us observe that we have the equalities dY(t) +Asgn(F(f )) df = dW(t) — 
psgn(F(f))dVi(f) + C7dW2(f) by virtue of (3.1), (3.3), (1.4), therefore 

dY{t) = sgn(y(f)) (pdVi(t)-Xdt)+adW 2 (t). 

But this is an equation of the form (8.7) with A(-) = 0, for which Proposition 8.1 (iii) 
is satisfied, and the Brownian motions Vi(-) and W2{ ) are independent. Thus, path- 
wise uniqueness holds for this equation and, since the equation admits a weak solu- 
tion, this solution is actually strong: 

fWcf lW ((), 0<r<~. 

Furthermore, we deduce 3^(0 ^ ^^''^(f) from the second equation in (3.3), 
therefore also 

$ V (t) C t V{ - W2) {t), 0<f<oo; 

it follows from the last two displayed inclusions and (4. 1 3) that we have 

d (W U W 2 )^ = d (Y.V) {t) Q d (VuW 2 )( t) Q d (WuW 2 )^ ( < f < oo, 

and this proves the complementability of Vi(-) by W2{ ) ■ As for the third claim of 
this paragraph, it is fairly straightforward from the strict inclusion in (4. 13) that V\ (•) 
cannot be complemented by V2G) m tne two-dimensional filtration F^' ^ 2 ' . 

Finally, let us turn to the Brownian motion V(-) = pV\ (■) + crV2(') in (3.4). We 
follow the same approach as for the Brownian motion Q (•) , only with the roles of 
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8 = 2 pa and y = p 2 — a 1 interchanged, and with W (•) replaced by £/(•) from (3.2). 
A bit more precisely, we observe 

7 f sgn(Y(s))dV{s) + 8U b (t) = y{pW x {t) - oW 2 {t)) + 8{aWi{t) - pW 2 (t)) = 
Jo 

= pWi(f) + aW 2 (t) = W(t) = Y(t)-y-X f sga{Y(s))ds, 

Jo 

or equivalently 

Y(t) = y + J^sgn(Y(s))d(rV(s)+Xs)+8U b {t), 0<f <°°. 

Since we are assuming 8 = 2pa > 0, it follows from Proposition 8.1 (iii) that 
pathwise uniqueness holds for the above equation (even when y = p 2 — a 2 = 0); 
and since the equation admits a weak solution, this solution is actually strong. To 
wit, the process Y(-) is adapted to F^ V ' U ' , and we arrive at the filtration identities 
F< W ^F(™') — p(w,v,£/ b ) _p(Wi,W 2 ) (recall Proposition 4.2). Thus U b (-) com- 
plements the Brownian motion V (•) in F^ 1,1 ^ . 
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